Thermal rectifying effect in two dimensional anharmonic lattices 
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We study thermal rectifying effect in two dimensional (2D) systems consisting of the Frenkel 
Kontorva (FK) lattice and the Fermi-Pasta-Ulam (FPU) lattice. It is found that the rectifying 
effect is related to the asymmetrical interface thermal resistance. The rectifying efficiency is typically 
about two orders of magnitude which is large enough to be observed in experiment. The dependence 
of rectifying efficiency on the temperature and temperature gradient is studied. The underlying 
mechanism is found to be the match and mismatch of the spectra of lattice vibration in two parts. 

PACS numbers: 67.40.Pm, 63.20.Ry, 66.70.+f, 44.10.+i 
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in recent years. After two decades analytic and nu- 
merical studies in ID model, much progress has been 
achieved. On the one hand, the study has enriched 
our understanding about the underlying physical mech- 
anism. On the other hand, the study has made it 
possible to seek the practical application of heat control 
and management. Indeed, in 2002, Terraneo et al^ 
proposed a thermal device which can rectify the heat 
current through it when reversing the temperature 
gradient. The model proposed by Terraneo et al is a 
ID anharmonic lattice consisting of three segments with 
the Morse on-site potential of different parameters. As 
the first attempt of controlling heat current, the ratio 
of the thermal current changes is less than two. More 
recently, Li et ali^ 5 . construct a thermal diode model in 
which two Frenkel-Kontoroval (FK) chains with different 
nonlinear strengths are connected by a harmonic spring. 
The most successful improvement of the model by Li 
et al2^ lies in three facts: First, the configuration is 
more simple, it consists of only two different segments; 
Second, the ratio of heat current from two different 
directions is increased drastically about 100 times; 
Third, the underlying mechanism of thermal rectifying 
effect is explained by illustrating the phonon bands 
of the particles in different segments. Following Li et 
al's work, we further improve the rectifying efficiency 
from 100 to 2000 by substituting the weak FK chain 
with a Fermi-Pasta-Ulam (FPU) chain 2 ^. In addition, 
we find that the rectifying effect (asymmetric heat 
flow) is closely related to asymmetric interface thermal 
resistance (also called Kapitza resistance). Moreover, a 
specific relationship between the ratio of heat current 



and the overlap of the phonon spectra is demonstrated 
numerically. The rectifying effect can also be achieved 
by modulating the periodicity of the on-site potential of 
the FK latticed. 

The above works demonstrate the possibility of con- 
trolling heat current by changing structures/parameters 
of anharmonic lattices. These might find potential ap- 
plication in energy saving material. However, almost all 
works so far arc focused on ID system of finite size. Ob- 
viously, much progress has been achieved, but the final 
purpose is to put these ideas to application. It is thus 
a nature step forward to seek effective thermal devices 
to control heat current in higher dimension. The open 
questions are: whether the heat control mechanism in ID 
is still valid for the high dimension(s) and whether the 
extra dimension(s) reduces the rectifying efficiency? The 
answer might not be trivial, as in higher dimension the 
lattice vibration includes not only the longitudinal one 
but also the transverse ones. The longitudinal modes 
will couple to the transverse modest, which might affect 
the rectifying efficiency. 

In this paper, we concentrate our study on a 2D recti- 
fier model. We will demonstrate with numerical evidence 
that a 2D rectifier can be built up and it works in a very 
wide temperature range. The 2D thermal rectifier shows 
similar behaviors with ID thermal rectifier under similar 
parameters regimes. 

The paper is organized as the follows. In Section II, 
we describe our model and numerical method used for 
the computer simulation. In Section III, we demonstrate 
and discuss the dependence of the thermal rectifying effi- 
ciency on the temperature and the temperature gradient. 
Section IV is devoted to the interface thermal resistance 
which is the key point to understand the thermal rectify- 
ing effect. In Section V, we give physical understanding 
of the rectifying effect in terms of the lattice vibration 
spectra (also called phonon band.) We conclude the pa- 
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FIG. 1: (Color Online) Configuration of a two dimensional 
(2D) thermal rectifier from the Frenkel Kotoroval and the 
Fermi-Pasta-Ulam lattices. The left part is a 2D FK lattice 
and the right one is a 2D Fermi-Pasta-Ulam lattice. The two 
parts are connected by harmonic springs with constant ki n t. 
The left and the right ends are put into contact with heat 
baths of temperature Tl and Tr, respectively. 

per by conclusions and discussions in Sec. VI. 

II. MODEL AND METHODOLOGY 

In our previous work2&, we construct a ID thermal 
diode model by connecting a FK lattice to a FPU lattice 
with a weak harmonic spring. We denote it as 1D-FK- 
FPU model. This model displays a very good rectifying 
effect. In this paper, we extend the ID- FK-FPU thermal 
rcctihcr model to two dimensional one. We denote it as 
a 2D-FK-FPU model. The configuration of the 2D-FK- 
FPU model is illustrated in FijfFJ The left part is a plane 
of harmonic oscillators on a substrate whose interaction 
is represented by a sinusoidal on-site potential. Here we 
plot the contour line of the 2D sinusoidal potential. For 
simplicity, we put one particle in each valley. The right 
part is an array of an-harmonic oscillators represented by 
the FPU model. The two parts are connected by weak 
harmonic springs. The Hamiltonian of the system can be 
written as: 

H = HpK + Hfpu + H int , (1) 

where Hfk, Hfpu, Hi nt is the Hamiltonian of the left 
part, the right part and the interface section, respectively. 
They are represented in ©, ©, respectively. 
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Hint = / ] kj n t Vfi(\rNFK,r,NFK+i,j 1 — fo) (4) 

where r"i,j-k,i = — Qk,l is the relative displace- 
ment between particles, labelled as and (k,l). 
Vh(x) = \x 2 , Vfpu(x) = \x 2 + jx 4 , Ufk{x,v) = 

-T2#F C0S (f x ) 

The mass of the particles is uniformly 1. l is dis- 
tance between nearest neighbors in equilibrium. The 
particle at the ith column and the jth row is labelled 
as The coordinate and momentum of this particle 

is <fij = (xijiUij) andptj = (p Xtj ,p ViJ ). In order to es- 
tablish a temperature gradient, the two ends of the planes 
are put into contact with two Nose-Hoover heat bathes^ 
with temperature Tp and Tp for the left end and the right 
end, respectively. Particles for i = l,j = 1,2, 3...Ny 
are coupled with heat bath of temperature Tp and par- 
ticles for i = Nx,j = 1,2, 3...Ny are coupled with heat 
bath of temperature Tr. We checked in ID case that 
the result does not depend on the particular heat bath 
realization. Fixed boundary condition is used along tem- 
perature gradient direction, denoted as X direction in 
this paper, namely, <f ,j = (0, j), qN x +i,j = (N x + 
Periodic boundary condition is applied in the Y direc- 
tion, namely, cfi.i — g^jVy+i (see Fig. 1). Under these 
boundary conditions, the system can be considered as a 
tube. The total number of particles is Nx x Ny- All re- 
sults given in this paper are obtained by averaging over 
N x 10 s (N > 2)stcps after a sufficient long transient time 
when a non-equilibrium stationary state is set up. The 
equations of motion of the particles are: 

1i<j = Pi, 3 ' 

\-& = 2 ^-^ ( 5 ) 

and the auxiliary variables £jj are described by the 
equations: 



= 1 3=1 

+ kFKVHidnj-iJ-ll - lo)) + UFK(Xi,j,Vi,j)) 



n F k Ny -2 here, T is the temperature of the heat bath (Tp or Tr), 

it V^c^V \ u \r (\^ I ; \ and Q is the parameter of coupling between the thermal 

f-f *r~t 2 (2) bath and the system. In this study, we set Q = 1 so 

that the response time of the thermostats, , is of the 

same order of the original time scale of the lattice. Our 

purpose is to study rectifying effect in two dimension and 

n x N y -2 the dependence of the rectifying efficiency on the system 

Hfpu — y^(— i temperature and the temperature gradient, so we don't 

i=NpK+1 j—i ^ ,„s attempt to search the optimum setting of parameters. 

. , T/ /i-. i j \ We choose the system parameters same as in ID FK- 

+ kFPUVFPU\ \n,j;i-i,j - «o) ^ OTT , , , f , , , j . in 

FPU model which has been tested as good one m ID 

+ k FPU VFPu(\ri,j-i,j-i\ - lo)) caS e, that is k FK = 1, A = 5 , l = 1, k FPU = 0.2 and 
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0.05. We set 



T L =T (1 + A), 
T R =T (1-A), 



(7) 



where —0.8 < A < 0.8 in our simulations. So we can 
simply denote T as the temperature added on the system 
and 2A = (Tr — Tl)/Tq as the normalized temperature 
difference of the system. 

The temperature used in our numerical simulation is 
dimensionless. It is connected with the true temperature 
T r through the following relation^, T r = — ^ — T, where 
m is the mass of the particle and b is the period of ex- 
ternal potential, ojq is the vibration frequency, ks is the 
Boltzman constant. For the typical values of atoms, we 
have T r ~ (10 2 — 10 3 ) 9 , which means that the room tem- 
perature corresponds to the dimensionless temperature 
T~ (0.1 - 1). 

The local temperature is defined as 



T %3 = 9 m (^,j) 



(8) 



where (...) stands for a temporal average. The local heat 
current J^j is defined as the energy transfer per unit time 
from the particle labelled as to the nearest particles 
along X direction. 
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d(V(\r l+hr ^\ - l ) + Vijn^.^l - lo)) 
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where k = kFK, kppu or k int , depending on the site along 
X direction. For a 2D-lattice, we treat only heat current 
flowing along the X-direction. We denote the current 
from the particles in ith section to the next section in 
the X-direction simply as Ji(Ji = Ylf=i The total 

current of the system is averaged over all sections, 
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(10) 



In our simulations, the fluctuations of temporal heat 
current through each section are all less than five per- 
cents. We use |J+/J_| as a rectifying efficiency, to de- 
scribe quantitatively the rectifying performance of the 
system. J + is the current when A is positive (heat flows 
from the FK part to the FPU part) and J_ is the current 
when A is negative (heat flows from the FPU part to the 
FK part). 

The model we used is a simple extension from one di- 
mension to two dimension, however the results in Fi^5] 
show that it truly demonstrates good rectifying effect 
on heat current. In Fig [2 we can see the visible differ- 
ence between J + and J_ in very wide temperature range. 
The difference varies from few times to several hundreds 
times. 
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FIG. 2: (Color Online) Heat current of the system with Ny = 
4 and Nx = 10 versus T . |A|=0.5 



III. DEPENDENCE OF RECTIFYING EFFECT 
ON TEMPERATURE AND TEMPERATURE 
DIFFERENCE 

In this section, we study the dependence of the system 
performance on the temperature change. Fig [3] shows 
the rectifying efficiency \J+/J-\ versus T . In Fij^a), 
we can see that there exists an optimum performance 
(OP) of the rectifying effect when changing temperature 
To. We can define the temperature for the optimum per- 
formance as Top- In Fi^3j), Fi^3}: and Fi^3ji we show 
the dependence of the ratio |J+/J_| on temperature un- 
der different conditions. We found that Top depends 
on the system settings along Y direction. In F1^3Jd, the 
number of particles along Y direction varies from 4 to 8 
and 16 while other settings are kept unchanged. We can 
see clearly that Top shifts to lower temperature when 
increasing N Y - The value of T OP is 0.04, 0.037, 0.0325 
for Ny = 4, 8 and 16, respectively. Top keeps the same 
value when we change Nx- This is shown in Fi^3h. In 
Fi^3]i, we change the periodic boundary condition in Y 
direction to free boundary condition. Top and the opti- 
mum performance change drastically. Top changes from 
0.037 to 0.025 and the ratio increases almost 100%. 

Quantity Wt in the figures is a parameter defined as 
the width of the effective temperature range T e over the 
half value of OP, while i? = Wt/Top is defined as the 
quality factor. T e , Wt, and d = Wt/Top are useful pa- 
rameters to estimate the temperature range in which the 
system has a good rectifying effect. In our investigation, 
T e broadens from the center [0.03-0.46] to high tempera- 
ture region or low temperature region. The typical value 
of Wt under different settings is around 0.018-0.023. i9 
is always larger than 0.5. The results suggest that the 
system is effective in very wide temperature range. 

In Fi^QJ we show the heat current versus the (normal- 
ized) temperature difference, A, for three different T . 
Full symbols represent J + and empty ones | J_ | . We can 
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FIG. 3: (Color Online) The rectifying efficiency, \J+/J-\, ver- 
sus To at different conditions. (a)The system with Ny = 4, 
N x = 10, and |A| = 0.5. Top = 0.04. (b) Compari- 
son for three different Ny with Nx = 40. Top = 0.04, 
W T = 0.022 and t? = 0.55 for the case of Ny = 4. When 
N Y = 8, Top = 0.037, W T = 0.022 and •& = 0.59. When 
N Y = 16, Top = 0.0325, Wt = 0.018 and ■d = 0.55.(c) 
Comparison for three different Nx with Ny = 8. They have 
the same Top = 0.037. When N x = 10, W T = 0.023 and 
= 0.62. When N x = 20, VFt = 0.020 and ■d = 0.54. 
When N x = 40, W T = 0.022 and = 0.59. (d) Com- 
parison for different boundary conditions along Y direction 
with Ny — 8, Nx = 20. Free boundary: The optimum 
performance (OP)OP = 298.5, T OP = 0.025, W T = 0.014, 
= 0.56; Periodic boundary: OP = 191.5, Top = 0.037, 
W T = 0.022, i? = 0.59. 
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FIG. 4: (Color Online) Heat current of the system with Ny = 
8, Nx = 20 versus temperature difference at three different 
To. Note that |J_| decreases as |A| increases in the interval 
0.1 < |A| < 0.6, this is the so-called "negative differential 
thermal resistance" , see context for more explanation. 



10 J 



♦10^ 



10 1 



10 u 



^-T =0.05 
T Q =0.045 
T Q =0.037 
slope=9.5 



0.0 



0.2 



0.4 



|A| 



0.6 



0.8 



FIG. 5: (Color Online) The ratio of heat current, |J+/J_| 
versus half normalized temperature difference |A[ at temper- 
ature T = 0.037, 0.045, 0.05. The dotted line has a slope of 
9.5. 



see that J+ increases with A monotonically and it always 
larger than \J-\, while J_ changes with A in different 
ways. One can find that there are three regions for \J—\. 
In the first region, |A| < 0.1, the increase of A leads to 
the increase of the heat current. However, in the second 
region, 0.1 < |A| < 0.6, the increase of |A| does not in- 
duce the increase of \J—\, instead it results in a decrease 
of \J-\. In the third region, A > 0.6, \J-\ is almost a 
constant independent of |A|. In this region, the \J- \ is so 
small that the system can be approximately considered 
as an insulator. The ratio |J_|_/J_| becomes larger and 
larger when |A| increases. That indicates that rectifying 



effect increases with increasing temperature difference. 

The strange behavior of |J_| observed in the second 
region, namely, the larger the temperature difference the 
smaller the heat current, is called negative differential 
thermal resistance. As we will demonstrate later that 
this is a typical phenomenon in nonlinear lattices. It 
can be understood from the match and mismatch of the 
vibrational spectra of the interface particles. 

Fi^5] shows \J+/J-\ versus |A|. It is found that 
|J_I_/J_| increases with |A| in an exponential way in the 
regime in which |A| is smaller than A2. Approximately, 

\J+/J-\ ocexp(c*|A|) + <5, (11) 



5 



1.5 



1.2 



0.9 



0.6 



=mm&=s= 



-■-Tq=0.05,A=0.5 

□ -Tq=0.05,A=-0.5 
• -Tq=0.07,A=0.5 
-o-Tq=0.07,A=-0.5 
a-Tq=0.09,A=0.5 
-a-Tq=0.09,A=-0.5 



12 



16 



20 



FIG. 6: (Color Online) T/T versus lattice site for T o =0.05, 
0.07, and 0.09. The solid symbols are for the cases of A = 0.5 
and the empty symbols are for the cases of A = —0.5. Ny = 
8,N X = 20. 



here c is about 9.5 under the particular parameter setting 
Ny = 8 and Nx = 20 with periodic boundary condition 
along Y direction. When we change Ny , Nx or boundary 
condition along Y direction, c changes slightly, but it is 
always around 10. The variation of |c — 10 1 is smaller 
than 0.5 in our investigation. S is related with Tq and A. 
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IV. INTERFACE THERMAL RESISTANCE 
(ITR) - KAPITZA RESISTANCE 

Thermal resistance between two different materials or 
between twin or twist boundaries of the same material 
has been extensively studied both experimentally and 
theoreticall y 30 ' 31 ^ 2 ' 33 . In fact, the existence of a ther- 
mal boundary resistance between a solid and superfluid 
helium was first detected by Kapitza 33 as early as 1940's. 
This boundary resistance is named Kapitza resistance af- 
ter him. Later, it is found that such a Kaptiza resistance 
exists at the interface between any pair of dissimilar ma- 
terials. KhalatnikovSi developed the acoustic mismatch 
model to explain the Kapitza resistance. Since then, con- 
tinuous efforts have been devoted to this problem. More 
information can be found in the review by Swartz and 
Pohl 3 ^. 

The Kapitza resistance is defined as 

R = AT/ J, (12) 

where J is the heat current and AT the temperature 
difference between two sides of the interface. 

In our system, the temperature drops are different 
when the temperature gradient of the system is reversed 
as shown in Fig. 6. Therefore, we use i?+ and i?_ to 
denote the interface resistance for the case of A > and 
A < 0, respectively. 



FIG. 7: (Color Online) (a) Interface thermal resistance R± 
versus Tq. |A| = 0.5. (b) The ratio R-/R+ versus T . N Y = 
4, N x = 10. 



In Fi^7] and Fi^SJ we show the dependence of IRT on 
temperature To and the normalized temperature differ- 
ence A. From the two figures we can see that, generally, 
R- (with larger temperature drop) is about two or three 
orders of magnitude larger than R + (with smaller tem- 
perature jump). Both R + and i?_ decrease with To until 
T) 0.2, then both become approximately constants. In 
Fifd 3 ! we show the ratio versus temperature Tq. 

It is clearly seen that there exists an optimal temperature 
value for the ratio R-/R+. In Fi^SJ we show the resis- 
tance versus temperature difference. We can see that i?„ 
monotonically increases with temperature difference un- 
til it reaches a maximum value, while R + monotonically 
decreases with increasing A. If we plot the ratio of i?_ 
over i? + versus temperature difference, we find that the 
relationship between them also obeys the exponential law 
like the ratio of heat current. 

Comparing Fi^T] and Fi^S] with Fi^3] and Fi^jSl we can 
find that the behaviors of the IRT and heat current 
through the system is very similar, both are asymmet- 
ric, both the ratio R-/R+ and J+/J_ have an optimum 
value under different temperature and obey exponential 
law when changing temperature difference. The asym- 
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FIG. 8: (Color Online) (a) R± versus the normalized temper- 
ature difference |A|. To = 0.037. (b) The ratio R-/R+ versus 
|A|. N Y = 8,N X = 20. 



metry of thermal resistance when reversing temperature 
gradient is the determinant factor for the rectifying effect 
on heat current of the system from the formula p2[) . 



V. PHYSICAL MECHANISM OF RECTIFYING 
EFFECT: AN ANALYSIS OF LATTICE 
VIBRATION SPECTRUM 

From above investigation, we know that the asymmet- 
ric interface thermal resistance determines the asymme- 
try of heat current when reversing temperature gradient 
on the system, but what cause the asymmetrical behav- 
ior of the fTR? In this section, we will answer this ques- 
tion from a fundamental point of view: lattice vibration 
spectrum. Lattice vibration is responsible to heat trans- 
port in our model. An effective way to get the spectrum 
of lattice thermal vibration is the discrete faster Fourier 
transform (DFFT) for the lattice velocity^. 

In our model, the vibration has two components, v x 
and v y . We can use the theorem of equipartition of en- 
ergy to simplify the numerical calculation. According to 
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FIG. 9: (Color Online) DFFT of v x ,v y . The shad- 
owed regions represent the integral of the power spectra. 
2S{P x {uj/2tv)/L) = 0.05 = (v 2 x ) « T,2S(P y (u/2w)/L) = 
0.05 = (Vy) « T . 



the equipartition theorem, the molecules in thermal equi- 
librium (here we have local thermal equilibrium) have 
the same average energy associated with each indepen- 
dent degree of freedom of their motion and that energy 
is ksT/2. For our system, we have m(v 2 )/2 = ksT, 
m{vl)/2 = m{v 2 y )/2 = k B T/2. 

In our calculation, m — 1, ks — 1. So we 



have {vl) = (v$) = T. If we do the DFFT of v x 
and v y separately, we should have, X^o* l^zWI 2 = 
A/ * Eiio _1 2|-Px(i)| 2 and E^" 1 \v y (t)\* = A/* 
X^=o _1 2|-Py(j)| 2 - Here P(f) is the Fourier transform 
of velocity v(t). From the equipartition theorem, one 
has Efjo 1 \v x (t)\ 2 /L = (vl) and Y.to M)? I L = (v 2 y ) 
From above analysis, we have A/^^q -1 2\P x (j)\ 2 / L = 

T, A/ J2f=o^ 2\P y (j)\ 2 /L = T. The power spectra of 
vibration obtained from DFFT of v x and v y agree with 
the above formula very well (see Fi^S])- The integral of 
power spectra is exactly equal to the temporal average 
of v 2 and v 2 individually. There is a slight difference be- 
tween (v x ) and (Vy). The difference might be caused by 
the number of sampled data or the different boundary 
conditions in X, Y direction. In the formula, L should 
be infinite. 

In our system, we find that the asymmetrical ITR and 
heat current are strongly related with the overlap of vi- 
bration spectra of the particles at the two sides of the 
interface. When the vibration spectra overlap with each 
other, the system behaves like a thermal conductor, while 
the system behaves like a thermal insulator when the vi- 
bration spectra are separated. 

Physically, whether an excitation of a given frequency 
can be transported through a mechanical system depends 
on whether the system has a corresponding eigenfre- 
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quency. If the frequency matches, the energy can easily 
go through the system, otherwise, the excitation will be 
reflected. In our system, the overlap of the two spectra 
means that there exists common vibrational frequency in 
two parts of the system. The excitation (here as phonon) 
of such common frequency can be transported from one 
part to the another. However, if the vibrational spec- 
tra of two parts are separated, then the excitation at 
any part cannot be transported to another part, because 
there exists no such corresponding frequency in another 
part. 

The change from overlap to separation is induced by 
the different temperature dependence of the vibration 
spectra of the two segments, which is a general feature 
of any anharmonic lattice. 

In Fi EfTUl we show the vibration spectra of the FK part 
and the FPU part under different temperatures. We can 
see that the vibration spectra of the FK part broaden 
from high frequency to low frequency when increasing 
the system temperature. This is because at low tem- 
perature, the atoms of the FK model is confined at the 
valley of the on-site potential, thus the atoms oscillate 
in very high frequency, however, when the temperature 
is increased, more and more low frequency modes can 
be excited. In the limiting case, when the temperature is 
larger enough that the kinetic energy of the atom is much 
larger than the on-site potential, then the FK model be- 
comes a chain of harmonic oscillators which has frequency 

On the contrary, the vibration spectra of the FPU part 
broaden from low frequency to high frequency. In fact, 
we have showr>2£ that the highest oscillation frequency of 
the FPU model depends on temperature, ujfpu ~ T 1 / 4 . 
Therefore, in some settings, the vibration spectra of the 
FK part and the FPU part will overlap with each other, 
while in other temperature settings, they will separate 
with each other. These are shown in FijTTI The compar- 
ison of vibration spectra of the FK part at two different 
temperature ranges and the FPU part at the full temper- 
ature range from 0.01 to 0.12 are shown in Fi jTTk and 
Fi jTTb separately. We can see that the vibration spectra 
of the FK part and the FPU part are matched with each 
other when the temperature of the FK part is from 0.05 
to 0.12 (see Fi efTTb ). The vibration spectra of the two 
parts are separated from each other when the tempera- 
ture of the FK part is from 0.01 to 0.03 (see FiflTb). This 
indicates that when the temperature of the FK part is be- 
low a certain value, which we call T c (0.04 ~ 0.045), the 
system will behave like a thermal insulator since the sep- 
arated vibration spectra of the interface particles make 
the heat conduction almost impossible. When the tem- 
perature of the FK part is above the critical point T c , 
the system will be a good thermal conductor since the 
matched vibration spectra allow the heat flow. Thus, if 
we adjust To and A appropriately to make Thigh > T c and 
Tiow < T c , the system will have a good rectifying effect. 
If T hlg h,Ti ow > T c or Thigh, Ti ow < T c , the rectifying 
effect is very poor. The above analysis is based on the 
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FIG. 10: (Color Online) Vibration spectra of the interface 
particles at different temperature, (a) The vibration spectra 
of the particles in the FK segment. (b) The vibration spectra 
of the particles in the FPU segment. 



vibration spectra of the system with Ny = 8, Nx = 20. 
It is consistent with the result obtained in Section III. 

Now we look back at Fi^3] The optimum point is 
at T = 0.037 for N Y = 8,N X = 20, corresponding 
T Mgh = 0.0505 > T c and T low = 0.0235 < T c . The effec- 
tive temperature range T e with |A| = 0.5 is from 0.0261 
to 0.0481. The low temperatures are all smaller than T c 
and all high temperatures are larger than T c . Both de- 
creasing To and increasing To in the outside region of T e 
lead the system to the two extreme cases Thigh ,Ti ow > T c 
or Thigh, Ti ow < T c with poor performance. From the 
above analysis, we can say that the different properties 
of heat current under different temperature and temper- 
ature difference are determined by the temperature de- 
pendence of the vibration spectra of the two segments. 

In terms of the vibration spectra of the particles in the 
interface, the complex behavior of |J_| in Fig 4 can be 
explained from the vibration spectrum. In particular, in 
the range from Ai to A2, a novel phenomenon- called 
the negative differential thermal resistance phenomenon 
is observed in Ref. 25 and fully discussed in Ref^I. In this 
particular temperature interval, a larger temperature dif- 
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FIG. 11: (Color Online) Comparison of the spectra of the FK 
part and the FPU part at different temperature region, (a) 
The spectra of the FPU part at T FPU = 0.01 - 0.09 (from 
bottom to top) and the FK part at T FK = 0.01 - 0.04 . (b) 
The spectra of the FPU part at T FPU = 0.01 - 0.09 and the 
FK part at T FK = 0.045 - 0.09 

ference can induce a smaller heat current. The negative 
differential thermal resistance can be understood from 
the overlap and separation of the vibration spectra of 
the interface particles. This phenomenon is valid for a 
wide range of the parameters. Moreover, in Rcf 37 , Li. 
et. al. show that it is this negative differential ther- 
mal resistance property that makes the thermal transis- 
tor possible. 

More importantly, we find a specific relationship be- 
tween the overlap of the vibration spectra of the two seg- 
ments and the ratio R-/R+ in the interface or the ratio 
| J + /J_| from two directions. We introduce the following 
quantity to describe overlap of the vibration spectra, 

q fpnf)p r x (f)df jpnf)p r x (f)df , . 

± JP?(f)dffP?(f)df T int T M 1 J 

S± corresponds to the case of A > and A < 0, re- 
spectively. In Fi dT^l we plot S + /S_ versus and 
|J + /J_|. A very good power law was found between the 
ratio of resistances or heat currents and the overlap of 




FIG. 12: (Color Online) The bullet are the ratios of heat 
current |J+/J_| versus S+/S- . Square are the ratios of the 
ITR R-/R+ versus S+/S- . The dotted line has a slope 1.35± 
0.03. 



the vibration spectra: \J + /J-\ ~ oc (S+/S-) 1 

The best fit for the ratio of current suggests the power 
law constant 7 = 1.35 ± 0.03. The two figures Fi efTTI and 
Fi Ell2| give us a very clear and quantitative picture about 
the dependence of the rectifying effect of the system on 
the vibration spectra. 

Since the rectifying effect sensitively depends on the 
overlap of the vibration spectra, we may find answers in 
Fi Jl3l for the behavior of the 2D system responding to the 
temperature changes at different conditions. We can see 
clearly that when we change the number of particles in Y 
direction, the temperature for the optimum performance 
will change, whereas when we change the number of par- 
ticles in X direction, the temperature for OP are kept at 
the same value. And the value of the temperature for the 
OP at different conditions found by the overlap are con- 
sistent with the value in Fi^3l These results suggest that 
the different vibration spectra of the two segments and 
the overlap between them are the determinant factors of 
the system complex behaviors. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper, we have studied the thermal rectify- 
ing effect in a 2D anharmonic lattice. The performance 
of the 2D thermal rectifier under different environment 
changes, such as system temperature and the tempera- 
ture difference on the two sides of the system, have been 
investigated systematically We find that there exits an 
optimum performance (OP) for a specific thermal rec- 
tifier at certain temperature range. The OP is affected 
by the boundary condition and the number of particles, 
Ny, along Y direction. The OP shifts to lower tem- 
perature when increasing Ny or changing the periodic 
boundary condition to the free boundary condition along 
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FIG. 13: (Color Online) The ratio ol the overlap of vibration 
spectra S+ / S- versus To at different conditions, (a) Compar- 
ison at three different Ny with Nx = 20. (b) Comparison at 
three different Nx with Ny = 8. 



mal rectifier is the temperature difference between the 
two ends. We find the rectifying efficiency increases ap- 
proximately as an exponential law in certain temperature 
range with the temperature difference. The rectifying ef- 
ficiency is mainly determined by the asymmetrical ITR. 
The study on the ITR shows the similar behavior with 
heat current. 

The behaviors of the ITR and heat current of the sys- 
tem are strongly correlated with vibration spectra of the 
particles beside the interface. The asymmetry behavior 
of ITR and heat current is induced by the different tem- 
perature dependence of the vibration spectra of the two 
parts beside the interface. We find the vibration spec- 
tra of the FK part broaden from high frequency to low 
frequency, conversely, the vibration spectra of the FPU 
part broaden from low frequency to high frequency as the 
temperature increases. The different temperature depen- 
dence of vibration spectra makes the system transition 
from a thermal conductor to an insulator possible by set- 
ting the system temperature and temperature difference 
properly. Moreover a specific relationship between the 
performance of the system and the convolution of the vi- 
bration spectra of the two parts is found numerically as 
power law. 

Our study on 2D thermal rectifier gives a very clear pic- 
ture about how the system responds to the environment 
changes. The results should be useful for further experi- 
mental investigation. The thermal diode constructed by 
a monolayer thin film or a tube like structure might have 
many practical applications. 
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